Chapter 27 



Irrationally winding 



I don't care for islands, especially very small ones. 
— D.H. Lawrence 

(R. Artuso and P. Cvitanovic) 

THIS CHAPTER Is conccmed with the mode locking problems for circle maps: 
besides its physical relevance it nicely illustrates the use of cycle expan- 
sions away from the dynamical setting, in the realm of renormalization 
theory at the transition to chaos. 

The physical significance of circle maps is connected with their ability to 
model the two-frequencies mode-locking route to chaos for dissipative systems. 
In the context of dissipative dynamical systems one of the most common and 
experimentally well explored routes to chaos is the two-frequency mode-locking 
route. Interaction of pairs of frequencies is of deep theoretical interest due to the 
generality of this phenomenon; as the energy input into a dissipative dynamical 
system (for example, a Couette flow) is increased, typically first one and then two 
of intrinsic modes of the system are excited. After two Hopf bifurcations (a fixed 
point with inward spiralling stability has become unstable and outward spirals to 
a limit cycle) a system lives on a two-torus. Such systems tend to mode-lock: 
the system adjusts its internal frequencies sUghtly so that they fall in step and 
minimize the internal dissipation. In such case the ratio of the two frequencies 
is a rational number. An irrational frequency ratio corresponds to a quasiperiodic 
motion - a curve that never quite repeats itself. If the mode-locked states overlap, 
chaos sets in. The likelihood that a mode-locking occurs depends on the strength 
of the coupling of the two frequencies. 

Our main concern in this chapter is to illustrate the "global" theory of circle 
maps, connected with universality properties of the whole irrational winding set. 
We shall see that critical global properties may be expressed via cycle expansions 
involving "local" renormalization critical exponents. The renormalization theory 
of critical circle maps demands rather tedious numerical computations, and our 
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Figure 27.1: Unperturbed circle map (k = O'm (27.1)) 
with golden mean rotation number. 
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intuition is much facilitated by approximating circle maps by number-theoretic 
models. The models that arise in this way are by no means mathematically triv- 
ial, they turn out to be related to number-theoretic abysses such as the Riemann 
conjecture, already in the context of the "trivial" models. 



27.1 Mode locking 

The simplest way of modeling a nonlinearly perturbed rotation on a circle is by 
1-dimensional circle maps x ^ x' = fix), restricted to the one dimensional torus, 
such as the sine map 

k 

Xn+i - fi^n) - x,i + Q sinilnxn) mod 1 . (27.1) 

In 

f(x) is assumed to be continuous, have a continuous first derivative, and a con- 
tinuous second derivative at the inflection point (where the second derivative van- 
ishes). For the generic, physically relevant case (the only one considered here) the 
inflection is cubic. Here k parametrizes the strength of the nonlinear interaction, 
and D. is the bare frequency. 

The state space of this map, the unit interval, can be thought of as the elemen- 
tary cell of the map 

Xn+i ^ fi^n) = Xn + ^ - ^ sin(27r.£„) . (27.2) 

where " is used in the same sense as in chapter 25. 
The winding number is defined as 



W{k,Q) - lim(x,j - XQ)ln. 



{113) 



and can be shown to be independent of the initial value xq. 
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For ^ = 0, the map is a simple rotation (the shift map) see figure 27.1 

x„+i = x,j + Q mod 1 , (27.4) 
and the rotation number is given by the parameter Q. 

W{k = 0, Q) - Q . 

For given values of Q and k the winding number can be either rational or ina- 
tional. For invertible maps and rational winding numbers W = P/Q the asymptotic 
iterates of the map converge to a unique attractor, a stable periodic orbit of period 

Q 

p{xi) = Xi + P, / = 0,l,2,---,e- 1 . 

This is a consequence of the independence of xq previously mentioned. There is 

also an unstable cycle, repelling the trajectory. For any rational winding number, 

there is a finite interval of values of O values for which the iterates of the circle 

map are attracted to the PjQ cycle. This interval is called the PjQ mode-locked exercise 27.1 

(or stability) interval, and its width is given by 

Ap/e = Q-^'"-'^ = n'^Q - a^fg . (27.5) 

where 0.'p^Q i^'pfg) denote the biggest (smallest) value of Q for which W{k, fl) = 
P/Q. Paiametrizing mode lockings by the exponent fi rather than the width A 
will be convenient for description of the distribution of the mode-locking widths, 
as the exponents // turn out to be of bounded variation. The stability of the P/Q 
cycle is 

Ap/e = ^= f'(^o)f(x,) ■ ■ ■ f\xQ.O 

For a stable cycle |Ap/q| lies between (the superstable value, the "center" of the 
stability interval) and 1 (the ^'pfg, ^pfg endpoints of (27.5)). For the shift map 
(27.4), the stability intervals are shrunk to points. As D. is varied from to 1, 
the iterates of a circle map either mode-lock, with the winding number given by 
a rational number P/Q e (0, 1), or do not mode-lock, in which case the winding 
number is irrational. A plot of the winding number W as a function of the shift 
parameter Q is a convenient visualization of the mode-locking structure of circle 
maps. It yields a monotonic "devil's staircase" of figure 27.2 whose self-similai^ 
structure we are to umavel. Circle maps with zero slope at the inflection point 
(see figure 27.3) 

/'(x,)-0, /"(x,)-0 
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Figure 27.2: The critical circle map (A; = 1 in (27.1)) 
devil's staircase [27.3]; the winding number W as func- 
tion of the parameter £1. 
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Figure 27.3: Critical circle map = 1 in (27.1)) with 
golden mean bare rotation number. 
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= 1, Xc = in (27.1)) are called critical: they delineate the borderline of chaos 
in this scenario. As the nonlineaiity pai^ameter k increases, the mode-locked 
intervals become wider, and for the critical circle maps = 1) they fill out the 
whole interval. A critical map has a superstable P/Q cycle for any rational P/Q, 
as the stability of any cycle that includes the inflection point equals zero. If the 
map is non-invertible (k > 1), it is called supercritical; the bifurcation structure of 
this regime is extremely rich and beyond the scope of this exposition. 

The physically relevant transition to chaos is connected with the critical case, 
however the apparently simple "free" shift map limit is quite instructive: in essence 
it involves the problem of ordering rationals embedded in the unit interval on a hi- 
eraixhical structure. From a physical point of view, the main problem is to identify 
a (number-theoretically) consistent hierarchy susceptible of experimental verifi- 
cation. We will now describe a few ways of organizing rationals along the unit 
interval: each has its own advantages as well as its drawbacks, when analyzed 
from both mathematical and physical perspective. 

27.1.1 Hierarchical partitions of the rationals 

Intuitively, the longer the cycle, the finer the tuning of the parameter Q required to 
attain it; given finite time and resolution, we expect to be able to resolve cycles up 
to some maximal length Q. This is the physical motivation for partitioning mode 
lockings into sets of cycle length up to Q. In number theory such sets of rationals 
are called Farey series. They are denoted by f^g and defined as follows. The 
Faixy series of order Q is the monotonically increasing sequence of all iiTcducible 
rationals between and 1 whose denominators do not exceed Q. Thus Pi/Qi 
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belongs to if < P, < Qi < Q and (PilQi) ^ 1. For example 
1 1 1 2 1 3 2 3 4 1 



_|1 112132341 
^ " 15' 4' 3' 5' 2' 5' 3' 4' 5' T 



A Farey series is characterized by the property that if Pi-i/Qi-i and Pi/Qi are 
consecutive terms of Tq, then 

PiQi-i - Pi~iQi = 1. 
The number of terms in the Farey series Fq is given by 



^ 3Q' 

3>(Q) = y m) = ^ + OiQlnQ). (27.6) 



Here the Euler function 4>iQ) is the number of integers not exceeding and rel- 
atively prime to Q. For example, 0(1) = 1, 0(2) = 1, 0(3) = 2, ...,0(12) = 
4,0(13) - 12,... 

From a number-theorist's point of view, the continued fraction partitioning of 
the unit interval is the most venerable organization of rationals, preferred already 
by Gauss. The continued fraction partitioning is obtained by ordering rationals 
corresponding to continued fractions of increasing length. If we turn this ordering 
into a way of covering the complementary set to mode-lockings in a circle map, 
then the first level is obtained by deleting A[i], Apj, • • • , A[f,j], • • • mode-lockings; 
their complement are the covering intervals £[,£2, ■ ■ ■ ,fai, ■ ■ ■ which contain all 
windings, rational and iiTational, whose continued fraction expansion starts with 
[a\, . . .] and is of length at least 2. The second level is obtained by deleting 
A[i,2], A[i,3], • • • , A[2,2], A[2,3], • • • , A[„,„,], • • • and so on. 

The nth level continued fraction partition <S„ = {^1^2 • • • a„] is defined as the 
monotonically increasing sequence of all rationals Pt/Qi between and 1 whose 
continued fraction expansion is of length n: 



Pi 1 
— - [ai,a2,---, an] 



Qi ^ " ^' 1 
ai H 



1 

02 + ■ ■ • — 
an 



The object of interest, the set of the irrational winding numbers, is in this partition- 
ing labeled by <Soo = {a\a2aj, ■ ■ •}, ajt e Z^, i.e., the set of winding numbers with 
infinite continued fraction expansions. The continued fraction labeling is particu- 
larly appealing in the present context because of the close connection of the Gauss 
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shift to the renormaUzation transformation R, discussed below. The Gauss map 



1 

Tix) = -- 

X 



1 

X 

, x = (27.7) 



([• • •] denotes the integer part) acts as a shift on the continued fraction representa- 
tion of numbers on the unit interval 



X = [ai,a2,a^, ■ . ■] — > T{x) - [a2,a^, ■ . .] . (27.8) 



into the "mother" interval Icnai...- 

However natural the continued fractions partitioning might seem to a number 
theorist, it is problematic in practice, as it requires measuring infinity of mode- 
lockings even at the first step of the partitioning. Thus numerical and experimental 
use of continued fraction partitioning requires at least some understanding of the 
asymptotics of mode-lockings with large continued fraction entries. 

The Farey tree partitioning is a systematic bisection of rationals: it is based 
on the observation that roughly halfways between any two lai^ge stability intervals 
(such as 1 /2 and 1/3) in the devil's staircase of figure 27.2 there is the next largest 
stability interval (such as 2/5). The winding number of this interval is given by the 
Farey mediant {P+P')I{Q+Q') of the parent mode-lockings PIQ and P' IQ'. This 
kind of cycle "gluing" is rather general and by no means restricted to circle maps; 
it can be attained whenever it is possible to arrange that the Qth iterate deviation 
caused by shifting a parameter from the correct value for the Q-cycle is exactly 
compensated by the Q'th iterate deviation from closing the 2'-cycle; in this way 
the two near cycles can be glued together into an exact cycle of length Q+Q' . The 
Farey tree is obtained by starting with the ends of the unit interval written as 0/1 
and 1/1, and then recursively bisecting intervals by means of Farey mediants. 

We define the «th Farey tree level T„ as the monotonically increasing sequence 
of those continued fractions {a\ ,a2,---, at] whose entries a,- > 1, i - 1,2, . . . ,k - 
1, ak >2, add up to Xf=i ai - n + 2. For example 

r2-|[4],[2,2],[l,l,2],[l,3]| =(^,^,^,^). (27.9) 

The number of terms in T„ is 2". Each rational in r„_i has two "daughters" in r„, 
given by 



[■■■,a] 

[■■■,a-\,2] [■■■,a + l] 
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Figure 27.4: Farey tree: alternating binary or- 
dered labeling of all Farey denominators on the nth 
Farey tree level. 




Iteration of this rule places all rationals on a binary tree, labeling each by a unique 
binary label, figure 27.4. 

The smallest and the largest denominator in r„ are respectively given by 

[n-2] = -L^, [1,1,...,1,2] = :^cxp" , (27.10) 
n-2 Fn+2 

where the Fibonacci numbers f „ are defined by f „+i = Fn+Fn-\', Fq = 0, Fi = 
1, and p is the golden mean ratio 



1 + V5 

p = — -— = 1.61803... (27.11) 



Note the enormous spread in the cycle lengths on the same level of the Farey tree: 
n < Q < p". The cycles whose length grows only as a power of the Farey tree level 
will cause strong non-hyperbolic effects in the evaluation of various averages. 

Having defined the partitioning schemes of interest here, we now briefly sum- 
marize the results of the circle-map renormalization theory. 



27.2 Local theory: "Golden mean" renormalization 

The way to pinpoint a point on the border of order is to recursively ad- 
just the parameters so that at the recurrence times t = n\,n2,nT,, - • • the trajectory 
passes through a region of contraction sufficiently strong to compensate for the 
accumulated expansion of the preceding «, steps, but not so strong as to force the 
trajectory into a stable attracting orbit. The renormalization operation R imple- 
ments this procedure by recursively magnifying the neighborhood of a point on 
the border in the dynamical space (by rescaling by a factor a), in the parameter 
space (by shifting the parameter origin onto the border and rescaling by a factor 6), 
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and by replacing the initial map / by the nth iterate /" restricted to the magnified 
neighborhood 



There are by now many examples of such renormalizations in which the new func- 
tion, framed in a smaller box, is a rescaling of the original function, i.e., the fix- 
point function of the renormalization operator R. The best known is the period 
doubling renormalization, with the recurrence times n, = 2'. The simplest circle 
map example is the golden mean renormalization, with recurrence times = f , 
given by the Fibonacci numbers (27.10). Intuitively, in this context a metric self- 
similarity arises because iterates of critical maps are themselves critical, i.e., they 
also have cubic inflection points with vanishing derivatives. 

The renormalization operator appropriate to circle maps acts as a generaliza- 
tion of the Gauss shift (27.38); it maps a circle map (represented as a pair of 
functions {g, /), of winding number [a, b,c,. . .] into a rescaled map of winding 
number {b,c, . . .]: 



Acting on a map with winding number {a, a,a,...], /?„ returns a map with the same 
winding number [a, a, . . .], so the fixed point of /?„ has a quadratic irrational wind- 
ing number W = [a, a, a, . . .]. This fixed point has a single expanding eigenvalue 
6a- Similarly, the renormalization transformation Ra^ . . .RajRcn = Raicn-a,, has 
a fixed point of winding number Wp = [ai,a2, . . . ,a,j,,ai,a2> ■ ■ ■]» with a single 
expanding eigenvalue 6 p. 

For short repeating blocks, 6 can be estimated numerically by comparing suc- 
cessive continued fraction approximants to W. Consider the P,-/Q,- rational ap- 
proximation to a quadratic irrational winding number Wp whose continued frac- 
tion expansion consists of r repeats of a block p. Let Q-r be the parameter for which 
the map (27.1) has a superstable cycle of rotation number P,-/Q,- = [p,p, ■ ■ ■ ,p]- 
The 6p can then be estimated by extrapolating from 



What this means is that the "devil's staircase" of figure 27.2 is self-similar under 
magnification by factor 6p around any quadratic in^ational Wp. 

The fundamental result of the renormalization theory (and the reason why all 
this is so interesting) is that the ratios of successive PrIQr mode-locked intervals 
converge to universal limits. The simplest example of (27.13) is the sequence of 
Fibonacci number continued fraction approximants to the golden mean winding 
number W = [1, 1, 1, ...] = ( VS - l)/2. 



fp{x) ^ Rfp{x) = af"{x/a) 




(27.12) 



{11. yi) 
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When global problems are considered, it is useful to have at least and idea on 
external scaling laws for mode-lockings. This is achieved, in a first analysis, by 
fixing the cycle length Q and describing the range of possible asymptotics. 

For a given cycle length Q, it is found that the narrowest interval shrinks with 
a power law 



Ai/e cx Q'^ (27.14) 



For fixed Q the widest interval is bounded hy P/Q = Fn-i/Fn, the nth con- 
tinued fraction approximant to the golden mean. The intuitive reason is that the 
golden mean winding sits as far as possible from any short cycle mode-locking. 

The golden mean interval shrinks with a universal exponent 

Ap/e « 2"^^' (27.15) 

where P = F„-i, Q = F„ and fii is related to the universal Shenker number 6\ 
(27. 13) and the golden mean (27. 1 1) by 



^1 = Ijli^ = 1.08218... (27.16) 
21np 



The closeness of /ii to 1 indicates that the golden mean approximant mode-lockings 
barely feel the fact that the map is critical (in the k^O limit this exponent is /i = 1). 

To summarize: for critical maps the spectrum of exponents arising from the 
circle maps renormalization theory is bounded from above by the harmonic scal- 
ing, and from below by the geometric golden-mean scaling: 

3/2 >//,„/„> 1.08218---. (27.17) 



27.3 Global theory: Thermodynamic averaging 

Consider the following average over mode-locking intervals (27.5): 



n(T) = 2 Z ^Ph- (27.18) 

2=1 {P\Q)=i 

The sum is over all irreducible rationals P/Q, P < Q, and Ap/Q is the width of the 
parameter interval for which the iterates of a critical circle map lock onto a cycle 
of length Q, with winding number P/Q. 
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The qualitative behavior of (27.18) is easy to pin down. For sufficiently neg- 
ative T, the sum is convergent; in particular, for r = -1, Q(-l) = 1, as for the 
critical circle maps the mode-lockings fill the entire Q. range [27.15]. However, 
as T increases, the contributions of the narrow (large Q) mode-locked intervals 
Ap/g get blown up to l/A^^g, and at some critical value of t the sum diverges. 
This occurs for t < 0, as £1(0) equals the number of all rationals and is clearly 
divergent. 

The sum (27.18) is infinite, but in practice the experimental or numerical 
mode-locked intervals are available only for small finite Q. Hence it is necessary 
to split up the sum into subsets <S„ = {/) of rational winding numbers Pi/Qi on 
the "level" n, and present the set of mode-lockings hierarchically, with resolution 
increasing with the level: 

Z„(T) = J] (27.19) 

ieS„ 

The original sum (27.18) can now be recovered as the z = 1 value of a "gener- 
ating" function n(z, r) = li„z"Z„(T). As z is anyway a formal parameter, and 
« is a rather arbitrary "level" in some ad hoc partitioning of rational numbers, 
we bravely introduce a still more general, P/Q weighted generating function for 
(27.18): 

OO 

a{q,T) = 2] 2 e~^^piQQ'^r,ipiQ _ (27.20) 
2=1 {P\Q)=l 

The sum (27.18) corresponds toq = 0. Exponents vp/q will reflect the importance 
we assign to the P/Q mode-locking, i.e., the measure used in the averaging over 
all mode-lockings. Three choices of of the vp/g hierarchy that we consider here 
con^espond respectively to the Farey series partitioning 

oo 

Q.{q,T)^Yj^^Qy' Q^""'"' ^ (27.21) 

2=1 (P\Q)=l 

the continued fraction partitioning 

oo 

a{q,T) = Y^e~''" 2 g2r^i„,....,„„] ^ (27.22) 

'1=1 ['3|,...,a„] 

and the Farey tree partitioning 

oo 2" 

a(q, T) = J] 2-^" ^i""' ' 2,VP,- e Tn . (27.23) 

k=n i= 1 

We remark that we ai^e investigating a set arising in the analysis of the parameter 
space of a dynamical system: there is no "natural measure" dictated by dynamics, 
and the choice of weights reflects only the choice of hierarchical presentation. 
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27.4 HausdorfF dimension of irrational windings 



A finite cover of the set irrational windings at the "nth level of resolution" is 
obtained by deleting the parameter values conesponding to the mode-lockings in 
the subset <S„; left behind is the set of complement covering intervals of widths 



.mm 



-a' 



.max 



(27.24) 



■Pr/Qr 



■Pi/Qi ■ 



Here ^'p"q (0'™^^^) are respectively the lower (upper) edges of the mode-locking 
intervals ^pjq^ (Api/Qi) bounding £j and / is a symbolic dynamics label, for ex- 
ample the entries of the continued fraction representation P/Q = [ai,a2, .-, an] of 
one of the boundary mode-lockings, / = aia2 • • • a„. provide a finite cover for 
the irrational winding set, so one may consider the sum 




The value of -t for which the n ^ oo limit of the sum (27.25) is finite is the 
Hausdorjf dimension Dh of the irrational winding set. Strictly speaking, this is 
the Hausdorff dimension only if the choice of covering intervals is optimal; 
otherwise it provides an upper bound to Dh- As by construction the ^, intervals 
cover the set of irrational winding with no slack, we expect that this limit yields 
the Hausdorff dimension. This is supported by all numerical evidence, but a proof 
that would satisfy mathematicians is lacking. 

The physically relevant statement is that for critical circle maps Dh = 0.870 . . . 
is a (global) universal number. exercise 27.2 

27.4.1 The HausdorfF dimension in terms of cycles 

Estimating the n — > oo limit of (27.25) from finite numbers of covering intervals 
£i is a rather unilluminating chore. Fortunately, there exist considerably more 
elegant ways of extracting Dh- We have noted that in the case of the "trivial" 
mode-locking problem (27.4), the covering intervals are generated by iterations 
of the Farey map (27.37) or the Gauss shift (27.38). The nth level sum (27.25) can 
be approximated by X", where 



£r(y,x)^dix-f-'(y))\f(y)\'' 



This amounts to approximating each cover width fj by \df" /dx\ evaluated on the 
ith interval. We are thus led to the following determinant 



det(l -z£^) 
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= nn(^-^""i^''i'/^y • ^^7.26) 

p k=0 



The sum (27.25) is dominated by the leading eigenvalue of Xr; the HausdorfF 
dimension condition Zni-Du) - 0(1) means that r - -Dh should be such that 
the leading eigenvalue is z - 1. The leading eigenvalue is determined by the 
k - part of (27.26); putting all these pieces together, we obtain a pretty formula 
relating the Hausdorff dimension to the prime cycles of the map f{x): 

^ = Y[(^-^I^^P^''") ■ ^27.27) 



For the Gauss shift (27.38) the stabilities of periodic cycles ai^e available analytical- 
ly, as roots of quadratic equations: For example, the Xa fixed points (quadratic 
irrationals with Xa - [a,a,a . . .] infinitely repeating continued fraction expansion) 
are given by 



x^=-"^^\ A.^-h^^] (27.28) 



and the Xab = [a, b, a, b, a,b,.. .] 2-cycles are given by 



2b 

2 _ (ab + 2+ \jab{ab + 4) ^ 

Kib - \XabXba) 



We happen to know beforehand that Dh - 1 (the irrationals take the full mea- 
sure on the unit interval, or, from another point of view, the Gauss map is not a 
repeller), so is the infinite product (27.27) merely a very convoluted way to com- 
pute the number 1? Possibly so, but once the meaning of (27.27) has been grasped, 
the corresponding formula for the critical circle maps follows immediately: 

0-]~[(l-l/|5/,|^") . (27.30) 



The importance of this formula relies on the fact that it expresses Dh in terms 
of universal quantities, thus providing a nice connection from local universal ex- 
ponents to global scaling quantities: actual computations using (27.30) are rather 
involved, as they require a heavy computational effort to extract Shenker's scaling 
5p for periodic continued fractions, and moreover dealing with an infinite alpha- 
bet requires control over tail summation if an accurate estimate is to be sought. In 
table 27. 1 we give a small selection of computed Shenker's scalings. 
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Table 27.1: Shenker's 6p for a few periodic continued fractions, from ref. [27.1]. 



p 


Sp 


[1111 ...] 


-2.833612 


[2222 ...] 


-6.7992410 


[3 3 3 3 ...] 


-13.760499 


[4 4 4 4...] 


-24.62160 


[5 5 5 5 ...] 


-40.38625 


[6666 ...] 


-62.140 


[12 12 ...] 


17.66549 


[13 13 ...] 


31.62973 


[14 14...] 


50.80988 


[15 15 ...] 


76.01299 


[2 3 2 3 ...] 


91.29055 



27.5 Thermodynamics of Farey tree: Farey model 



^Imr We end this chapter by giving an example of a number theoretical model 
motivated by the mode-locking phenomenology. We will consider it by means of 
the thermodynamic formalism of chapter L, by looking at the free energy. 

Consider the Fai^ey tree partition sum (27.23): the naiTowest mode-locked 
interval (27. 15) at the nth level of the Farey ti"ee partition sum (27.23) is the golden 
mean interval 

Af„_./f„ cx l^iT". (27.31) 

It shrinks exponentially, and for t positive and large it dominates q(T) and bounds 
dq{T)/dT: 

q'ma. = ^ = 1-502642... (27.32) 

However, for r large and negative, ^(t) is dominated by the interval (27. 14) which 
shrinks only harmonically, and q{T) approaches as 

^ ^ ^ ^ 0. (27.33) 
T ?iln2 

So for finite n, ^„(t) crosses the r axis at -t = D„, but in the n ^ oo limit, the 
q{T) function exhibits a phase transition; q{T) = for t < -Djj, but is a non-trivial 
function of t for -De < t. This non-analyticity is rather severe - to get a clearer 
picture, we illustrate it by a few number-theoretic models (the critical circle maps 
case is qualitatively the same). 

An approximation to the "trivial" Farey level thermodynamics is given by the 
"Farey model," in which the intervals {p/g are replaced by Q~^: 

2" 

Z«(T) = (27.34) 

!=1 
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Here Qi is the denominator of the ith Farey rational Pj/Q,. For example (see 
figure 27.4), 



Z2(l/2) = 4 + 5 + 5+4. 



By the annihilation property (27.38) of the Gauss shift on rationals, the nth Farey 
level sum Z„(-l) can be written as the integral 

Z„(-l) = j dx6(f"{x)) = Y,U\fa,...aM^ 

and in general 

Z„(t) = J dx£'^{0,x), 

with the sum restricted to the Farey level ai + . . . + ayt = " + 2. It is easily checked 
that /' „ (0) = (- 1 )*^ 2? 1 , so the Farey model sum is a partition generated by 
the Gauss map preimages of x = 0, i.e., by rationals, rather than by the quadratic 
irrationals as in (27.26). The sums are generated by the same transfer operator, so 
the eigenvalue spectrum should be the same as for the periodic orbit expansion, but 
in this variant of the finite level sums we can can evaluate q(T) exactly for t = k/2, 
k a nonnegative integer. First, one observes that Z„(0) = 2". It is also easy to check 
that Z„(l/2) = Zi Qi ^2- 3". More surprisingly, Z„(3/2) = Zi = 54 ■ T'K 
A few of these "sum rules" are listed in the table 27.2, they are consequence of 
the fact that the denominators on a given level are Farey sums of denominators on 
preceding levels. exercise 27.3 

A bound on Dh can be obtained by approximating (27.34) by 

Z„(t) = + 2V"^ (27.35) 



In this approximation we have replaced all {p/Q, except the widest interval {\/n, 
by the narrowest interval ^f„_,/F„ (see (27.15)). The crossover from the harmonic 
dominated to the golden mean dominated behavior occurs at the r value for which 
the two terms in (27.35) contribute equally: 



/ln«\ ^ In 2 

Dn = D + 0{ , D=-— = .72... (27.36) 

n ) Imp 



For negative t the sum (27.35) is the lower bound on the sum (27.25) , so D is 
a lower bound on Dh. 

From a general perspective the analysis of circle maps thermodynamics has 
revealed the fact that physically interesting dynamical systems often exhibit mix- 
tures of hyperbolic and marginal stabilities. In such systems there are orbits that 
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T/2 


Z„(r/2)/Z„_,(T/2) 





2 


1 


3 


2 


(5 + VT7)/2 


3 


7 


4 


(5 + VT7)/2 


5 


7 + 4V6 


6 


26.20249 . . . 



Table 27.2: Partition function sum rules for the Farey model. 

stay 'glued' arbitrarily close to stable regions for arbitrarily long times. This is a 
generic phenomenon for Hamiltonian systems, where elliptic islands of stability 
coexist with hyperbolic homoclinic tangles. Thus the considerations of chapter 24 
are important also in the analysis of renormalization at the onset of chaos. 



Resume 

The mode locking problem, and the quasiperiodic transition to chaos offer an 
opportunity to use cycle expansions on hierarchical structures in parameter space: 
this is not just an application of the conventional thermodynamic formalism, but 
offers a clue on how to extend universality theory from local scalings to global 
quantities. 



Commentary 

Remark 27.1 The physics of circle maps. Mode-loclcing phenomenology is re- 
viewed in ref. [27.6], a more theoretically oriented discussion is contained in ref. [27.3]. 
While representative of dissipative systems we may also consider circle maps as a crude 
approximation to Hamiltonian local dynamics: a typical island of stability in a Hamil- 
tonian 2-dimensional map is an infinite sequence of concentric KAM tori and chaotic 
regions. In the crudest approximation, the radius can here be treated as an external pa- 
rameter Q, and the angular motion can be modeled by a map periodic in the angular 
variable [27.9, 27.12]. By losing all of the 'island-within-island' structure of real sys- 
tems, circle map models slcirt the problems of determining the symbolic dynamics for a 
realistic Hamiltonian system, but they do retain some of the essential features of such 
systems, such as the golden mean renormalization [13.5, 27.9] and non-hyperbolicity in 
form of sequences of cycles accumulating toward the borders of stability. In particular, in 
such systems there are orbits that stay "glued" arbitrarily close to stable regions for arbi- 
trarily long times. As this is a generic phenomenon in physically interesting dynamical 
systems, such as the Hamiltonian systems with coexisting elliptic islands of stability and 
hyperbolic homoclinic tanglees, development of good computational techniques is here 
of utmost practical importance. 
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Remark 27.2 Critical mode-locking set The fact that mode-lockings completely fill 
the unit interval at the critical point has been proposed in refs. [27.3, 27.14]. The proof 
that the set of irrational windings is of zero Lebesgue measure in given in ref. [27.15]. 



Remark 27.3 Counting noise for Farey series. The number of rationals in the Farey 
series of order Q is (l>iQ), which is a highly irregular function of Q: incrementing Qhy I 
increases by anything from 2 to Q terms. We refer to this fact as the "Euler noise." 

The Euler noise poses a serious obstacle for numerical calculations with the Farey 
series partitionings; it blocks smooth extrapolations to Q ^ oo limits from finite Q data. 
While this in practice renders inaccurate most Farey-sequence partitioned averages, the 
finite Q HausdorfF dimension estimates exhibit (for reasons that we do not understand) 
surprising numerical stability, and the Farey series partitioning actually yields the best nu- 
merical value of the Hausdorff dimension (27.25) of any methods used so far; for example 
the computation in ref. [27 . 1 6] for critical sine map (27 . 1 ), based on 240 < Q < 250 Farey 
series partitions, yields Dh = .87012 ± .00001. The quoted error refers to the variation of 
Dh over this range of Q; as the computation is not asymptotic, such numerical stability 
can underestimate the actual error by a large factor. 



Remark 27.4 Farey tree presentation function. The Farey tree rationals can be 
generated by backward iterates of 1/2 by the Farey presentation function [27.17]: 

Mx) = x/(l-x) 0<x<l/2 

fi(x) = {l-x)/x l/2<x<l. l^/.J/; 



The Gauss shift (27.7) corresponds to replacing the binary Farey presentation function 
branch /o in (27.37) by an infinity of branches 



fa(x) = /i o/*"-"(x)= i-fl, _l_<x<i, 
fab-Ax) = fcO-ofboUx). " " (27.38) 



A rational x — {ai,a2, ■ ■ ■ ,ak\\s annihilated by the ^th iterate of the Gauss shift, /«,«,. (a:) - 
0. The above maps look innocent enough, but note that what is being partitioned is not 
the dynamical space, but the parameter space. The flow described by (27.37) and by its 
non-trivial circle-map generalizations will turn out to be a renormalization group flow 
in the function space of dynamical systems, not an ordinary flow in the state space of a 
particular dynamical system. 

The Farey tree has a variety of interesting symmetries (such as "flipping heads and 
tails" relations obtained by reversing the order of the continued-fraction entries) with as 
yet unexploited implications for the renormaUzation theory: some of these are discussed 
in ref. [27.13]. 

An alternative labeUng of Farey denominators has been introduced by Knauf [27.7] 
in context of number- theoretical modeling of ferromagnetic spin chains: it allows for a 
number of elegant manipulations in thermodynamic averages connected to the Farey tree 
hierarchy. 
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Remark 27.5 Circle map renormalization The idea underlying golden mean renor- 
malization goes back to Shenker [27.12]. A renormalization group procedure was for- 
mulated in refs. [27.8, 27. 10, 27. 11], where moreover the uniqueness of the relevant eigen- 
value is claimed. This statement has been confirmed by a computer-assisted proof [27. 18], 
and in the following we will always assume it. There are a number of experimental evi- 
dences for local universality, see refs. [27.19, 27.20]. 

On the other side of the scaling tale, the power law scaling for harmonic fractions 
(discussed in refs. [27.2, 27.13]) is derived by methods akin to those used in describing 
intermittency [27.24]: l/Q cycles accumulate toward the edge of 0/1 mode-locked inter- 
val, and as the successive mode-locked intervals l/Q, I HQ - 1) lie on a parabola, their 
differences are of order Q^^ . 

Remark 27.6 Farey series and tine Riemann liypotliesis The Farey series thermo- 
dynamics is of a number theoretical interest, because the Farey series provide uniform 
coverings of the unit interval with rationals, and because they are closely related to the 
deepest problems in number theory, such as the Riemann hypothesis [27.25, 27.26] . The 
distribution of the Farey series rationals across the unit interval is surprisingly uniform - 
indeed, so uniform that in the pre-computer days it has motivated a compilation of an en- 
tire handbook of Farey series [27.27] . A quantitative measure of the non-uniformity of the 
distribution of Farey rationals is given by displacements of Farey rationals for PilQi € Tq 
from uniform spacing: 

6i^— — , 1,2,---,(B(0 

O(0 a' 

The Riemann hypothesis states that the zeros of the Riemann zeta function lie on the 
s - 1/2 + IT line in the complex s plane, and would seem to have nothing to do with 
physicists' real mode-locking widths that we are interested in here. However, there is a 
real-line version of the Riemann hypothesis that lies very close to the mode-locking prob- 
lem. According to the theorem of Franel and Landau [27.28, 27.25, 27.26], the Riemann 
hypothesis is equivalent to the statement that 

Q,<Q 

for all e as 2 ~* The mode-lockings Ap/g contain the necessary information for 
constructing the partition of the unit interval into the £i covers, and therefore implicitly 
contain the 5, information. The implications of this for the circle-map scaling theory have 
not been worked out, and is not known whether some conjecture about the thermodynam- 
ics of irrational windings is equivalent to (or harder than) the Riemann hypothesis, but the 
danger lurks. 

Remark 27.7 Farey tree partitioning. The Farey tree partitioning was introduced in 
refs. [27.29, 27.30, 27.13] and its thermodynamics is discussed in detail in refs. [27.16, 
27.17]. The Farey tree hierarchy of rationals is rather new, and, as far as we are aware, 
not previously studied by number theorists. It is appealing both from the experimental 
and from the golden-mean renormalization point of view, but it has a serious drawback of 
lumping together mode-locking intervals of wildly different sizes on the same level of the 
Farey tree. 
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Remark 27.8 Local and global universality. Numerical evidences for global universal 
behavior have been presented in ref. [27.3]. The question was reexamined in ref. [27.16], 
where it was pointed out how a high-precision numerical estimate is in practice very hard 
to obtain. It is not at all clear whether this is the optimal global quantity to test but at 
least the Hausdorff dimension has the virtue of being independent of how one partitions 
mode-lockings and should thus be the same for the variety of thermodynamic averages in 
the literature. 

The formula (27.30), linking local to global behavior, was proposed in ref. [27.1]. 

The derivation of (27.30) relies only on the following aspects of the "hyperbolicity 
conjecture" of refs. [27.13, 27.21, 27.22, 27.23]: 

1. limits for Shenker 5's exist and are universal. This should follow from the renor- 
malization theory developed in refs. [27.8, 27.10, 27.11, 27.18], though a general 
proof is still lacking. 

2. 6p grow exponentially with «p, the length of the continued fraction block p. 

3. 6p for p = aiai- ■ .n with a large continued fraction entry « grows as a power 
of «. According to (27.14), lim„^oo(^p °^ n^ ■ In the calculation of ref. [27.1] the 
explicit values of the asymptotic exponents and prefactors were not used, only the 
assumption that the growth of 6p with « is not slower than a power of n. 

Remark 27.9 Farey model. The Farey model (27.33) has been proposed in ref. [27. 16]; 
though it might seem to have been pulled out of a hat, the Farey model is as sensible de- 
scription of the distribution of rationals as the periodic orbit expansion (27.26). 

Remark 27.10 Symbolic dynamics for Hamiltonian rotational orbits. The rotational 
codes of ref. [27.40] are closely related to those for maps with a natural angle variable, 
for example for circle maps [27.37, 27.39] and cat maps [27.41]. Ref. [27.40] also offers 
a systematic rule for obtaining the symbolic codes of "islands around islands" rotational 
orbits [27.43] . These correspond, for example, to orbits that rotate around orbits that rotate 
around the elliptic fixed point; thus they are defined by a sequence of rotation numbers. 

A different method for constructing symbolic codes for "islands around islands" was 
given in refs. [27.46, 27.44]; however in these cases the entire set of orbits in an island was 
assigned the same sequence and the motivation was to study the transport implications for 
chaotic orbits outside the islands [27.43, 27.45]. 



Exercises 

27.1. Mode-locked intervals. Check that when ^ the 
interval A/j/q have a non-zero width (look for instance 
at simple fractions, and consider k small). Show that for 



small k the width of Aq/i is an increasing function of k. 
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27.2. Bounds on Hausdorff dimension. By making use of 
the bounds (27.17) show that the Hausdorff dimension 
for critical mode lockings may be bounded by 

2/3 <Dh< .9240 . . . 

27.3. Farey model sum rules. Verify the sum rules re- 
ported in table 27.2. An elegant way to get a number of 
sum rules for the Farey model is by taking into account 
an lexical ordering introduced by Contucci and Knauf, 
seeref. [27.31]. 

27.4. Metric entropy of the Gauss shift. Check that 
the Lyapunov exponent of the Gauss map (27.7) is given 
by ;r^/61n2. This result has been claimed to be rele- 
vant in the discussion of "mixmaster" cosmologies, see 
ref. [27.33]. 

27.5. Refined expansions. Show that the above estimates 
can be refined as follows: 

F{z, 2) ~ ^(2) + (1 - z) logd - z) - (1 - z) 

and 

F(z, s) ~ ^(s) + r(i - .)(i - zy-' -s(s)(i - z) 



for s € (1,2) (S(s) being expressed by a converging 
sum). You may use either more detailed estimate for 
({s, a) (via Euler summation formula) or keep on sub- 
tracting leading contributions [27.34]. 

27.6. j„ and a^r- Look at the integration region and how it 
scales by plotting it for increasing values of n. 

27.7. Estimates of the Riemann zeta function. Try to 

approximate numerically the Riemann zeta function for 
i = 2, 4, 6 using different acceleration algorithms: check 
your results with refs. [27.35, 27.36]. 

27.8. Farey tree and continued fractions I. Consider the 
Farey tree presentation function / : [0, 1] i-» [0, 1], such 
that if / = [0, 1/2) and J = [1/2, 1], f\, = x/(l - x) and 
f\j = (1 - x)/x. Show that the corresponding induced 
map is the Gauss map g{x) - l/x - [l/x]. 

27.9. Farey tree and continued fraction II. (Lethal weapon 

II). Build the simplest piecewise linear approxima- 
tion to the Farey tree presentation function (hint: sub- 
stitute first the righmost, hyperbolic branch with a lin- 
ear one): consider then the spectral determinant of the 
induced map g, and calculate the first two eigenvalues 
besides the probability conservation one. Compare the 
results with the rigorous bound deduced in ref. [24. 17]. 
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